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Abstract 

We extend the recent study of the /c-body embedded Gaussian en- 
sembles by Benet et al. (Phys. Rev. Lett. 87 (2001) 101601-1 and 
Ann. Phys. 292 (2001) 67) and by Asaga et al. (|cond-mat/0107363 



and cond-mat / 1 736^ ) . We show that central results of these papers 



can be derived directly from the symmetry properties of both, the 
many-particle states and the random fc-body interaction. We offer 
new insight into the structure of the matrix of second moments of the 
embedded ensemble, and of the supersymmetry approach. We extend 
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the concept of the embedded ensemble and define it purely group- 
theoretically. 

PACS numbers: 05.45.+b, 05.30.Fk, 05.30.Jp, 03.65.Fd 
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1 Introduction 



Recently, novel results on the spectral properties of the /c-body embedded 
unitary and orthogonal ensembles of random matrices were obtained both 
for Fermions (Refs. [|[] and 0) and for Bosons (Refs. and [Q). We refer to 
these two pairs of papers as to ABRW. These results were based on explicit 
analytical calculations. The authors did mention that their work had its root 
in the underlying symmetry properties of the embedded ensembles. They 
did not, however, display these symmetry properties explicitly, nor did they 
derive their results with the help of such properties. 

In the present paper, we remedy this situation. We show that central 
results of ABRW follow from the symmetry properties of the many-particle 
states and of the random fc-body interaction. For simplicity, we confine 
ourselves to the fc-body embedded unitary ensemble EGUE(/;;) although our 
results can certainly be generalized to the orthogonal case. 

To define EGUE(A;), we consider m particles (Fermions or Bosons) in / 
degenerate single-particle states. The random fc-body interaction connects 
states in the Hilbert space of these m particles. It is useful to introduce, for 
arbitrary t, the dimension Nt of the Hilbert space of t particles. For Fermions, 
we have Nt = (^^ while for Bosons, Nt = (^~^^t~^^ • -^^^ simplicity, we write N 
for Nm- We show that for the embedded ensemble, three symmetry groups 
are relevant. These are (i) the group SU(Z) of unitary transformations of the 
degenerate / single-particle states; (ii) the group \J{Nk) of unitary transfor- 
mations of the A;-body interaction; (iii) the group U(A^) of unitary trans- 
formations of the Hilbert space. The group SU(/) governs the embedding. 
The group l]{Nk) is the symmetry group of EGl]E{k). The group U(iV) is 
the symmetry group of the Gaussian unitary ensemble (GUE), i.e., the free 
Gaussian unitary ensemble in dimensions. 

The matrix elements of the fc-body interaction in the space of m-particle 
states are Gaussian random variables with zero mean. Therefore, all in- 
formation necessary for calculating averages over EGUE(A;) is contained in 
the matrix of second moments of the ensemble. We show that the spec- 
tral representation (the "eigenvalue expansion" of ABRW) of the matrix of 
second moments constructed explicitly in ABRW can also be obtained from 
symmetry arguments. To this end, we decompose the random fc-body in- 
teraction into a sum of fc-particle basic interaction terms which transform 
according to the irreducible representations of SU(/). The latter are specified 
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by the bodyness (unitary rank) quantum number b introduced by Mon and 
French The eigenvalues appearing in the spectral representation of the 
matrix of second moments are uniquely defined by this quantum number, 
and the eigenvectors are given by SU(/) Clebsch-Gordan coefficients. 

We isolate the U(iV)-invariant part of the matrix of second moments. The 
result is unexpected as this part consists of two pieces. We find the expected 
term which has the form of the matrix of second moments of the free GUE, 
or, equivalently, of EGUE(m). However, we also find another term which 
has the form of the matrix of second moments of the degenerate Gaussian 
ensemble of multiples of the unit matrix, or, equivalently, of EGUE(O). We 
show that for /c -C m, the U(iV)-invariant part favors a Gaussian spectral 
shape and non-ergodic spectral properties. For 2k > m, on the other hand, 
that part favors a semicircular spectral shape, and ergodic behavior. 

We also investigate the invariance of the supersymmetric n-point gen- 
erating functions of the ensemble under the unitary transformations of the 
interactions, and we examine the consequences of that symmetry for the 
saddle-point solutions. Finally, we extend the concept of the embedded en- 
semble to cases where the embedding is governed by an arbitrary compact 
group different from SU(/). 

The paper is organized in the following way. Some basic relations for 
the representations of SU(/) in Fock space are presented in Section ^ The 
ensemble EGUE(/c), its symmetry and the SU(/) expansion of the embedded 
interactions are discussed in Section ^. The symmetry properties of the 
matrix of second moments are discussed in Section ^. In Section |^, we isolate 
the U(iV) -invariant part of the matrix of second moments and determine 
its relative importance. The properties of the supersymmetric generating 
functions are investigated in Section ^ Comments on the application of the 
ABRW approach to other ensembles and concluding remarks are contained in 
Section 1^. Some mathematical details, including the evaluation of eigenvalues 
of the matrix of second moments, are given in Appendices. We always treat 
Fermions and Bosons jointly. For the representation theory we follow the 
books by Wybourne and Lichtenberg 
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2 Fock Space and the Group SU(/) 



In this Section, we introduce the group SU(Z) and describe the action of group 
operations on composite operators. We use a running label i = 
for the degenerate single-particle states. We use second quantization and 
denote the single-particle creation and annihilation operators by af and a^, 
respectively. The vacuum state is written as | ) . We consider both Fermions 
and Bosons. Hilbert space is spanned by = N^n orthonormal m-particle 
states denoted by ImVm)- The index = 1,...,A^ is a running label. 
We emphasize that the basis states \mVm) are not necessarily simple states 
(Slater determinants for Fermions or product states for Bosons as used in 
ABRW): Linear combinations of such states are also admitted. No specific 
assumption on the choice of the basis will be made. The creation operators 
for the basis states are denoted by A^{mv„i)- 

An element u G SU(/) generates a unitary transformation of the single- 
particle states and a corresponding transformation 

t{u)4t+{u) = Yl '^i'i4 > u e su(0 , (1) 

i' 

of the creation operators af. The operators T{u) are generated by the trace- 
less linear forms of afaii. The transformation u induces a linear transforma- 
tion of the operators A^(mVm). These operators transform according to the 
irreducible representations D^''^{u), 

T{u)A^{mv^)T^{u) = ^ Dil^Ju)A-^{mv'J . (2) 

In the Fermionic (Bosonic) case, the matrices D^"^{u) are given, to within 
an equivalence transformation, by the totally antisymmetrized (totally sym- 
metrized, respectively) powers of u. In the sequel, we write the Young 
tableaux of the irreducible representations of SU(/) as h = (/ii, . . . , 
with hj denoting the number of columns of length j. The symbol 0-' stands 
for j zeros. In the case of Fermions, D^'^{u) belongs to the Young tableau 
fm = (0™^^10'^™^^); in the case of Bosons, it belongs to the Young tableau 
fm = (^TiO'^^). The annihilation operators A{mVm) transform according to 
the conjugate irreducible representation V^'^^u) = {D^'^^u))*, where fm = 
^Qi-m-i^Qm-i-^ for Fermions, and = (0'~^m) for Bosons. 
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In analogy with the famihar fractional-parentage technique of first quan- 
tization, we can expand the m-particle states \mvm ) into a sum of products 
of the fc-particle creation operators A~^{kvk) and the (m — A;)-particle states 
\m — kvm-k)- Explicitly, 

\^^rr.)= E A^kv,)\sV.)Cj:Zl^^ ■ (3) 

Because of the duality relation, Eq. ( p7D below, the particle rank k and the 
particle rank m — k are intimately linked. This is why we use here and in 
the sequel the label 

s = m — k. (4) 
The coefficients Cf"]""i are the SU(/) Clebsch-Gordan coefficients for the 

J k k J s s 

coupling of basic states of the represenations D^'' (u) and D^^ (n) to the basic 
states of the representation D-^™(-u) (in short, the Clebsch-Gordan coeffi- 
cients for the coupling {fkfs)fm)- In the fractional-parentage terminology, 
these SU(/) Clebsch-Gordan coefficients represent the fc-particle coefficients 
of fractional parentage of the m-particle states, C/^^/^^^ = [kvkSVs\]rnVm\- 
The origin of the combinatorial factor in Eq. (§) and further details are ex- 
plained in Appendix ^. 

The pair A^{kvk)A{kv'i?j of /c-particle interaction operators transforms ac- 
cording to the direct product of the irreducible SU(/) representations D^*'{u) 
and 'D^*'[u). With the help of an appropriate Clebsch-Gordan transforma- 
tion, this product can be reduced to a direct sum 

D^'^{u) X Vf'^iu) =J2D^'iu) (5) 

b 

of irreducible representations D^*'{u) of Young tableaux gi, . These are 
uniquely specified by the "bodyness" quantum number |^ h which assumes 
the values 6 = 0, . . . , fc. For Fermions, we have gi, = (O^^^IO'^^^^^IO^^^) while 
for Bosons, gt = {bO''~^b). (In keeping with ABRW, we confine ourselves for 
Fermions to the case of less than half filling, 2m < /.) The set of basic 
fc-particle interactions Bk{bwb) which transform according to the irreducible 
representations D^''{u) of SU(/), 

T{u)B,{hw,)T+{u) = T.K\n,M)Bkibw',) , (6) 
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is given by 

B.ibw,) = Y: A^kv,)Aikvi)Cf2f,., . (7) 

The SU(/) Clebsch-Gordan coefficients C^fl'^^f^y' accomphsh the reduction. 
The indices Wb label the individual states (rows) of the SU(/) representa- 
tion D^*'{u) of Young tableau gb- These running indices range from one to 
Mb = Nl — N^_i, the dimension of the representation D^^{u). The represen- 
tations of the Young tableaux gb are self-conjugate and integer. Therefore, 
we can assume that the D^^{u) are real and that the Bk{bwb) are Hermitean, 
B'^ihwb) = Bkihwb). No other assumptions about the choice of the irre- 
ducible representations D^'^iu) are made. 
Eq. implies the relations 



{Bk{bwb))k = SboyNk , {Bk{hwb)Bk{h'w'y) )k = Sbb'S^,^^'^^ . (8) 

The symbol {0)t denotes the trace of O in the space of t-particle states. 
According to the Wigner-Eckart theorem for the group SU(/), the matrix 
elements of the operators B^ipuib) with respect to the m-particle states 
\mVm) are products of the reduced matrix elements {m\\Bk{h)\\m) and of 
the Clebsch-Gordan coefficients C^J"^^ j , of SU(/) for the reduction of the 

direct product D^'^{u) x l)^"^{u) into the direct sum of the representations 

{mv^ \B,{bwb)\ mv'^) = {m\\B,mrn)C'fZlf..v'^ ■ (9) 

The definition @ of the reduced matrix elements generalizes the canonical 
definition of the reduced matrix elements of the irreducible tensor operators 
of SU(2) given by Racah to the present case. The dimension-dependent 

— 1/2 

factors are absorbed in the reduced elements. For further details on 

the Wigner-Eckart theorem for SU(Z), we refer the reader to the book by 
Lichtenberg |^ and the references given therein. 

The norm ( Bl{bwb) )m is independent of Wb and equal to the square of the 
reduced element. For later use we note that the norm of any properly nor- 
mahzed /c-particle operator Skip) of bodyness h has this same value: Making 
use of Eq. (|^), we find that every linear form Skip) of BkipWb) which satisfies 
the normalization condition {Sk{b)Sk{b))k = 1, obeys 

( S+{b)Sk{b) = ( m\\Bk{b) \\m)^ . (10) 
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The matrices (mvm \Bk{bwb) \ fnv'^ ) play an important role in the sequel. 
Therefore, it may be useful to give a physical interpretation of the concept 
of bodyness, and to establish the relation between these matrices and the 
corresponding quantities appearing in ABRW. As an example, we consider 
the simplest version of a pair of creation and annihilation operators: We take 
k = 1, i.e., we consider the pair alaj. For i ^ j, this operator has bodyness 
6 = 1: It moves one particle from the single-particle state j to a different 
single-particle state i. For i = j, on the other hand, the operator can be 
written as the sum of two terms with 6 = and 6=1, respectively. Indeed, 
we can decompose ajai into a traceless part a|aj — (1/0 I]j=i c^jC^j and the 
remainder, (1//) Sj=iaj%- By virtue of Eq. (||), the traceless part has bo- 
dyness 6 = 1, while the remainder has bodyness 6 = 0. By construction, 
the interaction terms of bodyness b = k cannot be simulated by an inter- 
action of lower particle rank and represent the genuine fc-body interaction 
which generically changes the state of not less than k particles. The example 
just considered suggests, and a more detailed consideration shows, that the 
matrices ( mVm \Bk{bwb) \ mv'^ ) bear a close relationship with the Hermitean 
matrices of ABRW. 'Indeed, the sets of indices (6^6) and (so) can be 
identified because the dimensions of both sets are the same and given by Mb- 
The differences between the matrices {mVm \Bk{bwb) \ mv'^) and the matri- 
ces C^" of ABRW are: (i) We admit any basis and drop the specialization 
to Slater determinants or product states employed by ABRW. (ii) The nor- 
malization condition imposed by ABRW on the matrices C^" differs from 
our Eq. (||). (iii) The matrices {mVm \Bk{bwb) \ mv'^) are introduced group- 
theoretically, while the matrices C*" were constructed explicitly. We return 
to this comparison in Section ^ below. 

3 The Embedded Ensemble 

After the preliminary steps of Section |^, we turn to a group-theoretical 
classification of the embedded ensemble EGUE(/c). This ensemble describes 
m identical particles distributed over / degenerate single-particle states which 
interact through a random fc-particle interaction of GUE type, 

W{k) = Y: W^^^.^{k)A-^{kv,)A{kvi) . (11) 
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The coefficients Wy^y'^{k) are independent Gaussian random variables 
^^fe<(^) = iWv'^vt,{k)y with moments 



W.,v'^ik) = , Wy^y>ik)W,^,>ik) = j^Kv',k.v'^ ■ (12) 

The overbar denotes the average over the ensemble. Obviously, we must have 
m>k. 

The analysis of EGUE(fc) is simplified when we express the interaction 
W{k) in terms of the operators Bkip Wh) introduced in Eq. (|^). For brevity 
we write n = bwb, with k = referring to the case 6 = 0. Making use of the 
orthonormality of the traces {Bk[K)Bk{n') )k (cf. Eq. (|D), we find 

W{k)=Y.Bt,{t,)Wu{>^) , Wk{^) = {W{k)B,{^))k . (13) 

K 

This expansion decomposes the interaction W{k) into parts of well-defined 
bodyness b with b = 0, . . . ,k. For / very large compared to m, the terms 
of highest bodyness b = k dominate and yield the main contribution to the 
average norm {W'^{k))m- 

The matrix representation of the ensemble in the m-particle space has 
the form 

Hy^y,Jk) = {mv^\Wik)\mv'^). (14) 

The dimension of the matrices H{k) is equal to the dimension N = Nm of 
the m-particle Hilbert space. We recall that we do not make any assumption 
about the basis except that the states \mvm.) are orthonormal. An explicit 
expression for the matrices H[k) can be found with the help of the fractional- 
parentage expansion Eq. (|^) of the m-particle states. With s = m — k, we 
find from Eq. ([TT|) (for more details we refer to Appendix |^) 

H..<ik) = H E (cfc- J* w^,,^{k) . (15) 

The dependence of the matrix elements Hy^^'^{k) of H{k) on the interaction 
matrix elements Wy^,yi^{k) is, thus, governed by the SU(/) Clebsch-Gordan 

coefficients C'/^^/^,,^ for the coupling {fkfs)fm- The expansion Eq. (|T^) in- 
troduces another parametrization of W{k) in terms of the coefficients Wk{f^)- 
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The new parameters Wk{n) are real and, by Eqs. ( ]I2| ) and (|T3|), are indepen- 
dent Gaussian random variables with moments 



A 



2 



Wki^) = , Wk{K)Wk{^') = —5,,, . (16) 

In the new parametrization, Hy^^i^lk) is given by 

Hvmv'S^) = J2i^'"m\Bk{i^)\mv'^)Wk{K) . (17) 

The average of any function F{H{k)) of the Hamiltonian H{k) over the 
ensemble is given by the integral 

F{H{k)) = J dfi{W{k))P{W{k))F{H{k)) , (18) 

where dn{W{k)) denotes the product of differentials of the matrix elements 
^vkv'^{k), and P(W{k)) the probability density 

{AT ^ / AT /2 

In the parametrization of Eq. (|1^, ( W^{k) )k = J2k Wki^)^ the measure 
dfi(W{k)) simplifies to the product of differentials of Wk{n)- 

The free Gaussian unitary ensemble in N dimensions (the GUE) is invari- 
ant under unitary transformations. More precisely, let U denote an arbitrary 
unitary matrix of dimension N, U & U(iV). Then, with H a GUE Hamil- 
tonian, H = UHU^ is also a member of the ensemble and appears in the 
ensemble with the same weight as the Hamiltonian H. EGUE(A;) does not 
have this symmetry. However, EGUE(fc) does possess the U{Nk) symme- 
try of the Gaussian unitary ensemble for k particles. Indeed, EGUE(/c) is 
invariant under the unitary transformation 

w.,Ak) = E U^,,,W,^,'{k){U+h., (20) 



of the interaction matrices Wvi_y'^{k) by any A^^ x A^^ unitary matrix U G 
l]{Nk). By this symmetry transformation, the Hamiltonian H{k) containing 
the coeficients W{k) is replaced by the Hamiltonian H{k) containing the coef- 
ficients W{k). The ensemble average remains invariant, F{H{k)) = F{H{k)). 
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For the special matrices U = D^''{u) which belong to the SU(/) repre- 
sentation Df''{u), the Hamiltonian H{k) takes the simple form H{k) = 
Df'^{u)H{k)[D^'^{u)Y' and is unitarily equivalent to H{k), a property which 
is lacking for a generic U G \J{Nk). For the parametrization Eq. ([I7D , the 
matrices H{k) are replaced by 

H.^v'Jk) = Y.{mv^\B,iK')\mvi)W,iK)Al,,iU) . (21) 

kk' 

Here A''{U) denotes the matrix representation of f/ G \J{Nk), carried by the 
matrices {kvk\Bk{n)\ kv'f. ) , 

This representation is unitary and real. It is the direct sum of two irreducible 
representations, the identity representation carried by the matrix of -Ba,.(0), 
and the (A^^l — l) -dimensional irreducible representation of the Young tableau 
(10^fe~^l) carried by the matrices of Bk^n) with 6 > 0. For U = Df''{u), the 
matrices A^{U) simplify to A^,^([/) = SwD^'^ More details on the 

b " 

matrices A^{U) are given in Appendix 0. The invariance of the ensemble 
under U{Nk) implies that whenever the integrals over Wkin) are performed 
and the ensemble averages are expressed in terms of the matrix elements 
of the operators Bk{K), the resulting formulae must be invariant under the 
replacement oi Bk{i^) by 

B,{k) = Y: Ai,,{U)B,{K') , U G U(iV,) . (23) 

k' 

As pointed out in ABRW, the properties of EGUE(A;) are closely related to 
those of the "dual" ensemble EGUE(s)=EGUE(m — k). The dual ensemble 
describes m particles interacting by an s-particle CUE interaction W{s). 
The Clebsch-Gordan coefficients C{"1^"} „ and C{™''7„ differ at most by a 

JkVkJsVs JsVsJkVk 

phase independent of the f 's. Therefore, the Hamiltonian -ff(s) of the dual 
ensemble can be written as 

ii^^^'M={^\Y.^C\:::,..^^ s = m-k. (24) 

We use the link established by duality between EGUE(A;) and EGUE(s) in 
Section ^. 
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4 The Matrix of Second Moments 



In ABRW, the matrix of second moments plays a central role for the analysis 
of EGUE(A;). In the present Section, we address this matrix from the point 
of view of group theory. 

The matrix elements Hy^^i^ {k) are Gaussian random variables. Therefore, 
all information necessary for the evaluation of averages over the ensemble 
EGUE(A;) is contained in the matrix of second moments, 

^v^v'^,v^v'^ik) = Hy^yi^{k)Hi,^i,i^{k) . (25) 

As shown by ABRW, central information on the properties of EGUE(fc) can 
be deduced if a "generalized eigenvalue expansion" for this matrix can be 
found. In particular, the application of Efetov's supersymmetric averaging 
technique to EGUE(A;) becomes possible. Such an expansion has the form 

a 

Here, the "eigenvalues" A° must be positive, and the "eigenvectors" C"^^, 
should be Hermitean. In ABRW, the expansion (^) was constructed explic- 
itly. 

We show now that this expansion is closely related to the bodyness ex- 
pansion of the Hamiltonian H{s) of the dual ensemble EGUE(s). Following 
ABRW, we note first that the second moments of the Hamiltonian matrix 
elements of the ensembles EGUE(A;) and EGUE(s) satisfy the "duality" re- 
lation 

NkH,^,,Jk)H,^,,Jk) = N,H,^,,Js)H,^,,Js) . (27) 

Here, the average on the left-hand side is over EGUE(fc), the average on the 
right-hand side is over EGUE(s). The duality relation can be verified by 
expressing the matrix elements of H{k) and H{s) in terms of the coefficients 
of fractional parentage C/^^7^„^ (Eqs. (|T5|) and (|2^) ), and taking the ensemble 
averages of the products of the interaction matrix elements Wy^y'^{k) and 
^vsv'X^)- The duality relation is robust with respect to a truncation of 
Hilbert space: It remains valid when some of the basis states \mvm.) are 
excluded. However, this relation is obviously violated when the random k- 
particle interactions are modified in such a way that their ensembles lose their 
unitary symmetry so that Eq. (|12]) is not valid any more. An important 
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example of such a modification is the restriction of the interaction to the 
terms of the highest possible bodyness [p!0[] . 

Another important observation relates to duality. The duality relation al- 
lows us to evaluate the second moments of EGUE(A;) in terms of the second 
moments of EGUE(s). This yields the moments Hy^yi^{k) Hy^y'^{k) as bilin- 
ear forms of matrix elements of operators Bs{k). From the U{Ns) symmetry 
of EGUE(s), the resulting expressions are invariant under the replacement 
of Bsin) by 

M^) = Y.K'.iU)B,iK') (28) 

k' 

for any U G \J{Ns). Here, A*(f/) denotes the transformation matrix de- 
fined analogously to Eq. (|22|) . The same "dual" U{Ns) invariance holds 
obviously for any ensemble average F{H{k)) calculated from the second mo- 
ments H^rnv'^{k)Hij^iii^{k) evaluated in this way. This concerns, in particular, 
the supersymmetric n-point generating functions of EGUE(A;). 

We turn to the matrix of second moments Hy^j}'^{s)Hi}^^'^{s) of the 
Hamiltonian H{s). We decompose the interaction W{s) into a sum of terms 
with well-defined bodyness, 

Hv^v'Js) = J2{mv^mK)\mv'^)WsiK) . (29) 

K 

We take the average of pairs of the random variables Ws{k) and express 
the matrix elements {mvm \Bs{ii) \ rnv'^^ ) with the help of the SU(/) Wigner- 
Eckart formula Eq. (|^). We use the duality relation Eq. (pTl) and get 

A.r.^'r.,^^.^Jk) = HZZik)m~Jk) = ± C'Z'h'^A\k)C'^Z'^ (30) 

bwf, 

where 

^^(^) = §-^{mmmm)^ , c^S, = c^r.l/.< • (31) 

The upper limit at the summation symbol indicates that the summation over 
b is restricted to 6 < s. Equation (|3l| ) constitutes the eigenvalue expansion of 
the matrix Ay^ff^^fj^y/^^k) , with A''{k) > the eigenvalues and = (C'^^^z ) 
the Hermitean eigenvectors. The eigenvalues A''(fc) are labelled only by the 
bodyness quantum number b. 
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We demonstrate the dual U{Ns) invariance of this spectral decomposition 
of A^mv'^,vmv'^i^) by writing Eq. (|30D as (in the sequel we suppress the upper 
limit s at the summation symbol) 

K 

Here V^lk) stands for 

V:^,,Jk) = ^iVA^Q^,^ = X^/N^,{mVra\B,iK)\mv'^) . (33) 

Under the U{Ns) symmetry transformations Eq. (pSD , the matrices V(A;) 
transform as 

V^ik) = Y: K'.iU)V^'ik) , U e ms) . (34) 

Since the matrix A^{U) is unitary and real, the sum over k appearing in 
Eq. ( p2D remains invariant under this transformation. This fact proves the 
invariance of A„„c;^,{)„<(/c). 

The expansion (|30|) gives the matrix of the second moments in terms 
of the SU(/) reduced matrix elements {m\\Bs{b)\\m) and in terms of the 
SU(/) Clebsch-Gordan coefficients C^/^'' 7 , . The reduced matrix elements 
of Bs{k) yield the eigenvalues and can be expressed in terms of the SU(/) 
recoupling coefficients {iifmfk)fsfk)fm\{fmifkfk)9b)fm)- However, as indi- 
cated by Eq. (0) , these matrix elements can be evaluated even more directly 
by constructing suitable interactions of the required bodyness and particle 
rank, and by calculating the corresponding trace. We present this calculation 
in Appendix y. In the Fermionic case this yields 

whereas in the Bosonic case we get 

According to Eq. (p7|), the eigenvalue expansions of the matrices of second 
moments of the two dual ensembles EGUE(fc) and EGUE(s) are related by 

E ^1.' A^(fc)C^L.' =NsY: C:^.A\s)C^^,, . (37) 
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When we compare these results with the eigenvalue expansion obtained 
in ABRW we see that, aside from normalization factors, the two expansions 
agree. Indeed, the eigenvalues found in ABRW agree with Eqs. ( P5[ ) and 
except for the factors X'^/N^. which are missing in ABRW. This difference is 
caused by the fact that in ABRW, the second moment of the interaction 
matrix elements is normalized to unity, whereas in Eq. (0) we have used 
the standard GUE normalization for k interacting particles. We conclude 
that the eigenvalues of ABRW are identical to the squares of the reduced 
matrix elements {m\\Bs{b)\\m), and that the eigenvalues C*^ are, aside from 
a factor ^/N due to the difference in normalization, given by the Clebsch- 
Gordan coefficients C^^^, . The factors A^^^ and Ng appearing on either side of 
the duality relations Eqs. ( pTj) and ( |37D are likewise absent in ABRW. Again, 
this is due to the difference in normalization of the random variables. 

5 Invariance under the Group U(A^) 

We have stressed in Section ^ that for the free GUE ensemble in dimensions 
(or, equivalently, for EGUE(m)), the matrix of second moments is invariant 
under the transformation of the Hamiltonian by any unitary N x N matrix 
U e U(A^), while for k < m, the matrix of the second moments of EGUE(/c) 
does not have this symmetry property. In the present section, we display the 
broken U(A^) symmetry of EGUE(/c) explicitly. 

To determine the U(A^)-invariant part A^^\k) of the matrix of second mo- 
ments A{k), we consider a transformation U G U(A^) with H{k) — >■ H{k) = 
UH{k)U~^. Under this transformation, the components C^^^, of the eigen- 
vector transform like the components of an U(A^) tensor with one covari- 
ant index and one contravariant index. Starting from Eq. (|30|), separating the 
U(A^)-invariant part of C^^j^i^C^^^,^, and working out the sum over k with 
the help of the duality relation Eq. (0), we find that the U(A^)-invariant 
part A^^\k) has the form 



NA%k) 

m - 1 



( 




1 - S{k) 




(38) 
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where 



1 NsA^i-s) 1 {H{k)y 



A»W^-(/^W>., (39) 

The U(A^)-invariant part of A{k) is thus specified by the square A°(fc) of the 
average spectral width and by the ratio S{k). This ratio was introduced in 
ABRW as a measure of the ergodicity of the spectral centroids. We remark 
in parentheses that the converse is also true: The average spectral width 



A^{k) and the ratio S{k) are governed entirely by A'^^^k). There are no 
contributions either to A^{k) or to S{k) arising from the U(iV)-non-invariant 
part of A{k). 

The result displayed in Eq. (^) is somewhat unexpected. Indeed, the 
U(A^)-invariant matrix A^^\k) is the sum of two terms, A^^\k) = A'^{k) + 
A^{k). The first, 

AG (k) = —AG(k)6^ 5c ^, (40) 

has the expected form: The Kronecker symbols carry the same indices as 
in the matrix of second moments of the free Gaussian unitary ensemble of 
N X N matrices or, equivalently, in EGUE(m). The average spectral width 
corresponding to the CUE described by this term is given by 

A^{k)=A'{k)N'^-^^. (41) 

The form of the second term, 

<^'^,^^.'Jk) = A^ik)5.„^^,J,^,^^, (42) 

is unexpected. Indeed, the indices on the Kronecker deltas occur in the same 
way as in the matrix of second moments of the degenerate Gaussian ensemble 
of multiples of the N x N unit matrix or, equivalently, in EGUE(O). The 
average spectral width corresponding to the degenerate Gaussian ensemble 
described by this term is given by 

A^(fc)=A°(fc) ^'^f_j^^ (43) 
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The ensemble EGUE(A;) thus differs from the free GUE ensemble not only 
by the presence of the U(A^)-non-invariant part but also in the structure 
of the U(iV)-invariant part A^^^k). This part contains, in addition to the 
expected term yl*^(/c), also the term A^{k). The relative importance of the 
terms A'^i^k) and A^{k) is governed by the parameter S{k). For fixed m and 
I, this parameter decreases with increasing k from the value unity at = to 
the value A^'^ at = m, with S = N''^ at A; = m/2. At A; = 0, the GUE term 
A'^{k) is absent since A^(0) = 0. With increasing k, the width increases, 
and for k > m/2 and large N, becomes dominant, A'^(A;) = A^{k) up to terms 
of order 0(A^~^). For finite ^ 1 and small k, the U(iV)-invariant part 
A^^\k) favors a Gaussian spectral shape and non-ergodic behavior, while for 
k > m/2, the favored spectral shape becomes semicircular and the spectral 
fluctuations, ergodic. 

The U(A^) symmetry of EGUE(/c) is obviously violated most strongly 
when k is close to m/2. A convenient measure for the degree of symmetry 
breaking is given by the ratio P{k) of the norm of the U(A^)-invariant part 
and of the norm of A{k). From Eqs. (|8|) and (|0|) we find that P{k) 
is given by 

^ Tr([AW(A:)]+AW(A:)) ^ 1 - 2S{k) + N'S\k) 
^ ' TT{A+{k)A{k)) {N^ -l){R{k)/2) ■ ^ ^ 

Here R{k) denotes the coefficient 



first introduced in ABRW as the measure of ergodicity of the spectral widths. 
The ratio P{k) is symmetric about the point k = m/2 and attains its mini- 
mum there. 

The results of the present Section cast new light on the most surprising 
result of ABRW. There it was found that for k <^ m, the spectral fluctua- 
tions of EGUE(/c) differ markedly from Wigner-Dyson form and tend towards 
Poissonian behavior. While the present analysis does not yield direct infor- 
mation on this question, the fact that for /c -C m the embedded ensemble 
EGUE(fc) is dominated by the U(A^)-invariant term A^[k), lends additional 
plausibility to this result. 
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6 Symmetry Properties of the 

Supersymmetric Generating Functions 

We show how symmetry arguments are involved in the use of the supersym- 
metry aproach introduced by Efetov ||^ and developed by Verbaarschot, Wei- 
denmiiller and Zirnbauer[]ll|. We use the same conventions as in Ref. 



We address first the one-point function and its saddle-point solution and 
then turn briefly to the generating functions of higher order. 

The ensemble average of the one-point function Z{E, H) of energy E is 
given by the graded integral 



Z{E, H)= j d/i(V') exp|#+ X 1 + J) 7/; - X H)ipf^ 



(46) 



over the field ip of 2N components ipavm- The components with a = are 
ordinary complex variables, the components with a = 1 anticommute. The 
matrices which appear in the exponent are written as direct products of the 
2x2 matrices acting on the indices a and the N x N matrices acting on 
the indices Vm- The complex energy = E + irj contains the infinitesimal 
term = 0+ introduced to assure convergence, and J is the source matrix 
Jav,^,a'v'^ = Saa'{-T^^ ^vmv'^- FoT uotatioual simphcity, we suppress the 
index k throughout and write H{k) = H. 

All information about the ensemble is contained in the matrix of second 
moments A. We first consider a toy model. We omit the U(A^) -non-invariant 
part of the matrix A. Then, A is given by the U(iV)-invariant part A^^^ = 
A'^ + A^ defined in Section |^. Making use of the transformation 

exp|-V(^+^)2| = (^^y^ /dtexp|-^(t^ + 2ztA°^+^)| , 

(47) 

we find that in this case, the function Z{E, H) simplifies to the convolution 

(^Y' Jdte.p{-^t^)WW^TJ^) (48) 



of the Gaussian of width vA° with the one-point function Z'^{E,A'~') of 
the free GUE in dimensions with the average spectral width a/A^. The 
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average level density p(°)(£') corresponding to Z^'^\E) has the width -\/A° 
with A° = A*^ + A^. As for the kurtosis of the average level density, we 
follow the convention of ABRW: Their quantity Q^^^ equals unity in case the 
spectral shape is Gaussian, and zero if the spectrum has semicircular shape. 
We find Q(°) = 1 - (A^/A°)2. Thus, g(°) decreases from the value one at 
A; = to zero aXk = m. For large A^, the average level density is given by the 
convolution of a Gaussian of width VA° with the semicircular distribution 
of width v^A^. In general, the actual behavior of Z[E, H) is much more 
complex than that of this toy model. However, the toy model may perhaps 
be not far from the truth for k <^m and for k close to m. In both cases the 
non-invariant part of A is comparatively small. 

We now address the one-point function in its full generality. With the 
help of the eigenvalue expansion of the matrix of second moments A and 
the Hubbard-Stratonovich transformation, the quartic term appearing in 
the exponent of Eq. ( P6| ) can be simplified and the integrals over can be 
performed. This yields 



(49) 



Z{E,H) = |dMa)exp|-^Ar^ 

In (^E+l X 1 + J-^(T^ X y^^^l 



where d/i((T) denotes the measure for integration over the graded 2x2 ma- 
trices a'^ introduced in the Hubbard-Stratonovich transformation. We use 

= V NA^C^. The angular brackets denote graded traces. Under a dual 
unitary transformation U G U{Ns), the matrices are replaced by the ma- 
trices introduced in Eq. (p^. Since A''(f/) are real unitary matrices, the 
transformation coefficients can be absorbed in the new integration variables 
a" = X^k' ^L'(^)^'^'' "^ith = EkIct")^- The integral thus remains 

invariant as required. 

The same symmetry consideration limits the form of the saddle-point so- 
lution. Indeed, with a^^ the saddle-point solution, the dual U(A^s) invariance 
implies that the saddle-point approximation to the generating function 



Z,^{E,H) = exp |- ^ In (^E+1 x 1 + J - J] <p x ^ | 



(50) 



must be invariant under the replacement of V by V, Eq. (|34D, for any 
U G \J{Ns). This condition can be fulfilled only when a^^ = for all 6 > 0. 
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The saddle-point solution found in ABRW has just this form, with 



(To = 



E 



2VA° \ 



E 



.2VA0 



• 1 . 



(51) 



We stress that this saddle-point solution exists only by virtue of the fact that 
the basis states \'mVm ) represent a complete set of states of the irreducible 
representation D-^"'{u). This can be seen in the following way. As discussed 
in ABRW, the saddle -point equation has the form 



E 



El X 1 



(52) 



Substitution shows that cr^p = ^mCq is a solution only if the traces of the 
matrices vanish for all 6 > 0. These traces differ only by a factor from 
the traces {Bm{n))m- Moreover, only scalar quantities can have a nonzero 
trace in the complete representation space of an irreducible representation. 
Hence this condition is met. However, the exclusion of some of the states 
ImVm) from the Hilbert space of the ensemble would make this argument 
invalid. The supersymmetry approach could then still be applied, but a 
saddle point would not exist generically. 

For the generating functions of higher order, the symmetry properties of 
the ensemble manifest themselves in the same way. As shown in detail in 
ABRW, the two-point function can be written as the graded integral 

Z(E + oj,E-oj,H) = Jd^i{a)expl^~Nj2{{a^)') 

- (\n (^(E ■ 1 + u+L) X 1 + J-X^fj'^ X y'^j^l (53) 



over the 4x4 graded matrices a'^ 



ap,a'p' J 



The doubling of dimen- 



sions of the graded matrices reflects the fact that we are dealing with two 
propagators, the retarded propagator of energy E + uj {p = 1), and the ad- 
vanced propagator of energy E — uj {p = 2). The diagonal matrix L with 
Lap,a'p' — Saa'Spp'{—y'^^ distinguishes between the two cases. The matrix 
J contains the source parameters. As discussed thoroughly in ABRW, the 
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integral is dominated by the saddle-point manifold 



^SD = ^boT VqT 



E 



sp 



2VA0 



1-i 



\ 



E 



2VA° 



■ L 



(54) 



Here, T = (Tap^a'p') denote the 4x4 graded matrices belonging to the coset 
space U(l, 1/2)/ [U(l/1) x U(l/1)]. The integral as well as the saddle-point 
approximation are manifestly invariant with respect to the dual symmetry 
transformations U G U(A^s), Eq. (|3^). Saddle-point manifolds of the same 
structure dominate also the graded integrals which represent the generating 
functions of higher order. 



7 Discussion and Conclusions 

In this paper, we have shown that some of the central results of ABRW can 
be deduced with the help of symmetry considerations. This fact explains 
why analytical progress in understanding the properties of EGUE(A;) has 
been possible at all. In addition, our work offers a deeper insight into the 
structure of EGl]E{k). We mention, in particular, the decomposition in 
Section ^ of the matrix of second moments into a part which is invariant 
under U G U(A^) and another which is not. 

We mention the following generalizations of our work. 

(i) The theory developed in Sections ^ - § for EGUE(/i;) is based on 
the SU(/) expansion formula Eq. (]T5|). We emphasize that the approach 



developed in ABRW and in the present paper is not restricted to this case 
but is much more general. Indeed, it applies likewise to ensembles defined by 
the same formula but with the SU(/) Clebsch-Gordan coefficients replaced 
by the Clebsch-Gordan coefficients of another symmetry group G. The will 
now discuss this case in order to illuminate the algebraic structure of the 
approach used in the present paper. 

We introduce a generalized ensemble by considering two independent sys- 
tems labelled j = k and j = s whose basic states \fjVj) transform according 
to the irreducible representations D^^{g) of dimension A^^ of the group G. 
We assume for definiteness that the group G, (7 G G, is a compact simple 
Lie group; the irreducible representations of G will be labelled by their high- 
est weights /, their basis states by the running indices v. We assume that 
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a non-trivial interaction W{k) of GUE type occurs only in system k, with 
^^'fc<(^) ~ {fkVk\W{k)\fkv'i.) Gaussian and distributed according to Eq. (JT^). 
The embedding of this interaction into a space of different dimension is ac- 
complished by projecting the product states \ fk'i^k)\fsVs) onto the subspace 
of states which transform according to the irreducible representation D^"^{g) 
contained in the direct product D^'^{g) x D^^{g). We denote the associated 
projection operator by I{fm)- The ensemble is then defined in terms of the 
Hamiltonian 

H{k) = IiU)Wik)IiU) . (55) 

Eq. ( |55D comprises the essence of the group-theoretical extension of the idea 
of an embedded ensemble. This extension is independent of the existence of 
Fermions and Bosons. It relies only on group-theoretical concepts. 
Using the standard composition formula 

\rmlmVm) = E Ih^k )\fsVs )C):itZ ' (56) 

we find for the matrix elements of H{k) 

Hr^^^y^.'^ (k) = E iCf:Jrf:zy (^) • (57) 

Here C'/^f/J^ denotes the Clebsch-Gordan coefficient of the group G for 
the coupling {fkfs)fm- The multiplicity index distinguishes different 
D^"^{g) of the same highest weight fm which may appear in the reduction of 
D^''{g) X Df"{g). The dimension of the matrices H{k) is N = fiN^, where 
denotes the dimension of D^'^{g) and /i the multiplicity with which this 
representation appears in the reduction. By construction, the ensemble is 
invariant under unitary transformations U G U(A^fc) of dimension A^^ of the 
interaction matrices W^^yi^ik). 

In analogy to Eq. (p^, we can introduce the dual ensemble 

H{s) = I{fm)W{s)I{U , (58) 

with the random GUE interaction W{s) acting now on the second system s. 
Eqs. ( pTf) and (|12D and analogous equations valid for the Hamiltonian H[s) 
and the interaction W{s) imply for the matrices of second moments A{k) 
and A{s) of the two dual ensembles the duality relation (see Eq. (^)) 
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To derive the eigenvalue expansion of the matrix of second moments, 
we expand the interaction W{s) in terms of a complete set of A'^^ basic 
interactions Bs{k) = B^{n) normalized according to {Bs{k,)Bs{k,'))s = 
with Bs{0) denoting again a normalized multiple of the unit operator. Here, 
{Os)s denotes the trace of Os in the Hilbert space of states | fgVs) of the 
second system. This yields 



W{s) =Y,Bs{k)Ws{k) , Ws{K)Ws{^i') = —5^^, . (60) 

Using this expansion and proceeding as in Section ^, we find the eigenvalue 
expansion 

Basic properties of the ensemble can be derived from this formula. The eigen- 
value expansion can eventually be simplified by adapting the choice of the 
operators Bs{n) to the group chain l]{Ns) D D^". We show in Appendix 
that the traces {Bs{k) of the operators Bs{k) in the Hilbert space of the 
composite system vanish for n ^ 0. This makes it possible to use the su- 
persymmetry approach in a meaningful way. The graded integrals which 
represent the n-point functions of the ensemble are dominated by the saddle 
points and/or saddle-point manifolds analogous to those discussed in the 
preceding section. 

(ii) We may consider EGUE(A;) as a member of a family of ensembles 
obtained by a modification of the GUE, the free Gaussian unitary ensemble 
in N dimensions. Eq. shows that the GUE in N dimensions, GUE = 
EGUE(m), may be written as 

^vmv'^ = J2i rnv„,\B^{bwb)\mv'^ )W{bwb) , (62) 

with the sum over b running over all b = 0, . . . , m and all corresponding 
Wb, and with W{bwb) Gaussian distributed with mean value zero and second 
moment ^ 

W{bwb)W{b'w'b,) = . (63) 
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When we restrict the sum over b to b not larger than k and renormahze the 
matrix elements by the 6-dependent factors K{b) = \J N/ Nk{m\\Bk{b)\\m) , 
the free Gaussian ensemble turns into EGUE(/c). Other restrictions of the 
sum over b and other renormalizations K{b) of the matrix elements will create 
other "modified CUE ensembles" of N x N matrices 

= T.'i^^nr\Bmibwb)\mv'^)K{b)Wibw,) . (64) 

bw,, 

The prime at the summation symbol indicates the restriction of the sum 
over b. By construction, the resulting modified ensembles are all invariant 
under the SU(/) transformation Ti. Df"^{u)T-l[D^"^{u)]'^ . On expressing 
the matrix elements {■mVm\Bm{bwi))\'mv'^) in terms of SU(/) Clebsch-Gordan 
coefficients and performing the appropriate SU(Z) recoupling, we find for the 
matrix of moments of the modified ensemble the eigenvalue expansion 

Av^v'^,V^v'm = ^VmV'JI^VmV'r^ = J2 ^v'Jl,^'' Cl'^Jl^ ■ (65) 

bWb 

Here, denote the SU(/) Clebsch-Gordan coefficients Eq. (pID, and A'' 

the eigenvalues 

iV - V ^^^b 
b 

The angular bracket expression stands for the SU(/) recoupling coefficient 
for the indicated change of the coupling scheme, with denoting the highest 
weight of the identity representation. For A'' positive, the properties of the 
modified ensemble thus can be studied by applying Efetov's supersymmetric 
averaging technique in the same way as for EGUE(A;). 

The cases mentioned under (i) and (ii) are very instructive. They show 
that our group-theoretical approach to EGUE(fc) is very general. Morerover, 
the results obtained for EGUE(A;) are typical for a wide class of embedded 
ensembles defined in terms of a general compact group G. 
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A The composition formula for A^(mVm) 

The single-particle creation operators af transform according to the SU(/) 
representation D^^^ (u) . The reduction of the direct product D-^j-^ [u) x D-^^ (u) 
contains the representation D^^{u) just once. Therefore, the creation opera- 
tors A^{mVm) can be written as 

A-'i^v^) = 4=^ E < • • • <c(:Z > (67) 

where C(^^]^ stands for the coupling coefficient 

^n-.-im ^flilfli2 f2V2fli3 ' ' fm-lVm-lflim ' \ / 

V2V3...Vm-l 

This coefficient is given in terms of products of the SU(/) Clebsch-Gordan 
coefficients Cjf_i^^_^j-^i^ for the coupling with j = 2,...,m. In 

the Fermionic (Bosonic) case, the coefficients C(^^]^ are totally antisymmet- 
ric (symmetric) functions of The factor (m!)~^/^ takes care of 
the proper normalization {\A{mVm)A^{mVm)\) = {iTiVml^Vm) = 1: The 
total number of contractions contributing to the norm is equal to ml. This 
statement applies equally to Fermions and to Bosons. 

Indicating the coupling scheme explicitly, we can rewrite Eq. (^) in the 
symbolic form 

A+{mv^) = -^{{a+a+y^ . . . a+Yz ■ (69) 
Modifying the coupling scheme we find 

VkVs 

(70) 

where C f^v^f^vs SU(/) Clebsch-Gordan coefficients for the coupling 

{fkfs)fm- Substituting this relation into Eq. (p9| ) and introducing the cre- 
ation operators A'^{kvk) and A'^{sVs), we obtain 

AHmvJ = J^T. A^kv,)A^sv,) Cfc- . (71) 

This is the composition formula which underlies Eq. (^). 



25 



Taking into account that, due to the anticommutativity (commutativity) 
of the Fermionic (Bosonic) operators af, the products A'^{kvk)A'^{sVs) can 
be coupled only to operators belonging to the irreducible tableaux fm, i.e., 

T.^^ikv,)A^isv,)Cj:,^^^,^ = for f^U, (72) 

and making use of orthonormality of Clebsch-Gordan coefficients, we find 
that the composition formula can be inverted. 



A+{kvk)A+{sVs) = 




/^From this equation it follows that the matrix elements of the creation op- 
erators A~^{kvk) have the form 

{mv^\A+{kv,)\sv,) = (^^ ' {Cj-llJ* . (74) 
Substituting this equation in 

{mVm\W{k)\mv'^) = ^ {mVm\A^{kvk)\sVs) {sVs\A{kv'^)\mv'^)W,,^^'^{k) 

(75) 

yields the formula Eq. ([T5|). 



B The matrices A^(?7) 



Multiplying Eq. (^21) by {kv'j^\Bk{K")\kvk) and summing it over Vk and 
with the help of the orthonormality relation (H), we arrive at the following 
explicit expression for A^/^(f/) : 

A''.(^)= E ^..^.(A;^,|5,(K)|A:5^)(t/+),^.^(K|5fe('^')I^^A')- (76) 

We substitute {kvk\Bk{S^)\kv'^,) = 6yi^yi^{Nk)~^^'^ , use Eq. (fTBD, and find that 

K.iU) = AUU) = 5.0 . (77) 
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The representation A^{U) thus has a fully reduced form: Under the unitary 
transformations U € U(A^fc), the matrix of -Bfc(O) remains invariant whereas 
the matrices of the Bj.{n) with 6 > transform like the components of an 
irreducible U(A'^fc) tensor of dimension A^| — 1 and of Young tableau (10^'=~'^1). 
For U = D^^{u), from Eq. it follows that 



= Y.K'^S^){kvu\B,{bw',)\kv',) . (78) 



For U = Dfk{u), the matrices A*^([/) simplify to A^,^(?7) = 5w D^^, ^^^{u) . 

C The reduced matrix elements {m\\Bs{h)\\m) 

To calculate the squares of the reduced matrix elements, we apply the ap- 
proach by Mon and French P] (cf. also Refs. [0, |13|). At the end of Section ^ 
it was shown that the value of {m\\B s{h)\\m)'^ can be obtained from the norm 
( S'^{h)Ss{h) )m of a suitable s-particle operator Ss(h) of bodyness h. Any such 
operator can be written as the product 



5^(6) =( \S,{h) (79) 




of the SU(/) invariant polynomial of the particle number operator given by 
the first factor and of a 6-particle operator of bodyness b denoted by Ship). 
The form of the polynomial is uniquely determined by the particle ranks 
of Sgip) and Sij{h)\ Since Sgip) is an s-particle operator, the polynomial 
has to vanish when acting on the states of t = 6, . . . , s — 1 particles. This 
reduces the choice of Sgip) to the choice of Si,{h). The 6-particle operators 
of bodyness h change the state of h particles and cannot be simulated by 
the operators of lower particle rank. For Fermions, a suitable choice then is 
Ship) = ^af . . . a^ttb+i . . .a2b while for Bosons, we use Ship) = ^{afY{a2Y /b\ . 
Here, ^ denotes a normalization factor. This yields 
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for Fermions and Bosons, respectively. We substitute this result in 



— 

K s - b , 



{Snb)S,ib))r 



(81) 



and fix the value of ^ from the normalization condition {S^{b)Ss{b) )s = 1. 
This yields 



{m\\BM\mf = {Stib)Ssib)), 
^rri — b^^ — m + k — b^ 

fm — b\ fl + m + b — l\ 



(Fermions) 
(Bosons) , 



(82) 



in keeping with Eqs. (P5|) and 



D The traces {Bs{i^)) 



m 



Using the composition formula (150) we find that 



E 

rmVmVkVsV'g 



(83) 



We express the Clebsch-Gordan coefficients with the help of the symmetry 
relation (cf. Refs. |, |ll) 



fkVkfsVs 



IN, 



(84) 



Here, fk denotes the highest weight of the irreducible representation con- 
jugate to D-f''{g); Vk the basis state of this representation; C'^fjfjf^y^ the 
Clebsch-Gordan coefficient for the coupling {fkfm,)fs] the multiplicity in- 



dex for this coupling; Uy^ll and 1/^^^"^"'^ the unitary matrices 



ifk) 



•Ok-^k 



^^k^fkVkfkVk ' 



JjUkfsfni) 



rsfsVs r^rmfmVm /TjUk) Y 

MM ^ ^fkVkfruVm fkVkfsVs^'^VkVk) 



c 



(85) 



(86) 
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In Eq. (|85|), (00) denotes the (fv) labels of the identity representation. Sub- 
stituting Eq. ( p4[) and performing the sum yields 

N N 
{Bs{k) U = y{Bs{k) )s = 5«o-^ . (87) 



The traces are thus nonzero only for k = 0. 
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